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ABSTRACT 

The discovery of a relationship between supermassive black hole (SMBH) mass and 
spiral arm pitch angle (P) is evidence that SMBHs are tied to the overall secular evo- 
lution of a galaxy. The discovery of SMBHs in late-type galaxies with little or no bulge 
suggests that an underlying correlation between the dark matter halo concentration 
and SMBH mass (Mbh) exists, rather than between the bulge mass and Mbh- In 
this paper we measure P using a two-dimensional fast fourier transform and estimate 
the bar pattern speeds of 40 barred spiral galaxies from the Carnegie- Irvine Galaxy 
Survey. The pattern speeds were derived by estimating the gravitational potentials of 
our galaxies from ifg-band images and using them to produce dynamical simulation 
models. The pattern speeds allow us to identify those galaxies with low central dark 
halo densities, or fast rotating bars, while P provides an estimate of Mbh- We find 
that a wide range of Mbh exists in galaxies with low central dark matter halo den- 
sities, which appears to support other theoretical results. We also find that galaxies 
with low central dark halo densities appear to follow more predictable trends in P 
versus de Vaucouleurs morphological type (T) and bar strength versus T than barred 
galaxies in general. The empirical relationship between Mbh and total gravitational 
mass of a galaxy (Mtot) allows us to predict the minimum Mtot that will be observa- 
tionally measured of our fast bar galaxies. These predictions will be investigated in a 
subsequent paper. 

Key words: galaxies: kinematics and dynamics; galaxies: haloes; galaxies: spiral; 
galaxies: structure; galaxies: fundamental parameters (pattern speed) 



1 INTRODUCTION 

The connection between the overall morphology and dynam- 
ics of galaxies and the centrally located supermassive black 
holes (SMBHs) they harbor provides a fundamental con- 
straint in the study of galaxy formation and evolution. One 
well-known example of this link is seen in the correlation 
between SMBH mass [Mbh) and the large-scale stellar ve- 
locity dispersion [a) of the bulge/spheroid encompassing it 
(Gebhardt et al. 2000; Ferrarese & Merritt 2000). Likewise, 
the discovery of a relationship between Mbh and spiral arm 
pitch angle (P; Seigar et al. 2008) suggests that SMBHs are 
tied to the overall secular evolution of a galaxy. Seigar et al. 
(2005, 2006) have already shown that there is a strong cor- 
relation between P and rotation curve shear rate (S), where 
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small values of P correspond to large S (see their Figure 3) 
and a large central mass concentration. More recent stud- 
ies have further indicated that SMBHs are linked to secular 
processes (e.g. Cisternas et al. 2011, Orban de Xivry et al. 
2011). 

The stellar velocity dispersion of the bulge/spheroid is 
commonly linked to Mbh (Gebhardt et al. 2000; Ferrarese 
& Merritt 2000). For example, M33 is a bulgeless late-type 
galaxy with the lowest nuclear BH mass ever measured with 
an upper mass limit of 1500 Mq (Gebhardt et al. 2001). 
The discovery of SMBHs in the centers of late-type galaxies 
with little or no bulge (Satyapal et al. 2007, 2008) suggests 
that the size, or mass, of SMBHs may be tied to the dark 
matter halo virial mass, or concentration (see, e.g., van den 
Bosch et al. 2007 for a relation). This hints at the cause of 
the AIbh — P connection and the reason why M33 fits in this 
relation. Comparisons of the corrected central velocity dis- 
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persion with the maximum rotation velocity of disk galaxies 
also reveal a good correlation between Mbh and the mass of 
the dark matter halo (Ferrarese 2002). Additionally, bulge 
size and dark halo mass show a direct correlation in general, 
but the relationship exhibits a lot of scatter (see, e.g.. Ho 
2007). Recently, Booth & Schaye (2010) have shown through 
simulations that Mbh appears to be linked to halo binding 
energy rather than the halo mass, with the central halo con- 
centration playing a role in the relationship's scatter. Other 
recent papers have also argued both for (Seigar 2011; Volon- 
teri et al. 2011) and against (Kormendy & Bender 2011) a 
coupling of SMBH masses and the properties of the dark 
matter halo. 

The secular evolution of a barred spiral galaxy is chiefly 
influenced by the rotation rate of the non-axisymmetric com- 
ponent (e.g., Lynden-Bell & Kalnajs 1972; Kalnajs 1991; Ko- 
rmendy & Kennicutt 2004). Simulation models have shown 
that after a bar initially forms, the pattern speed (f2p) re- 
mains fast when embedded in a dark matter halo with a low 
central density (Debattista & Sellwood 2000). This is due 
to reduced dynamical friction between the two components 
that would otherwise cause a rapid decrease in Q,p. Q,p also 
determines the locations of important resonance regions in 
galaxies with a single perturbation mode, such as a bar or 
grand design spiral. This allows the rotation of a bar to be 
described by the ratio TZ = RcR/Rbar, where Rcr is the 
corotation, or 1:1, resonance (CR) radius and Rtar is the 
bar semimajor axis length. When 1.0 < TZ < 1.4 a bar is 
deemed to be a fast rotator, while TZ > lA implies a slow 
rotating bar. Theoretical arguments by Contopoulos (1980) 
state that self-consistent bars must have TZ > 1.0. 

About 50% of all spiral galaxies display either a full, 
broken, or partial ring shaped pattern in their light distri- 
bution (Buta & Combes 1996). Although rings are observed 
in some non-barred galaxies, they are typically associated 
with barred galaxies due to the ease at which they form 
through bar-driven gravity torques (see Grouchy et al. 2010 
and references therein). Theoretical arguments suggest that 
inner rings in barred spiral galaxies are features of the in- 
ner ultraharmonic, or 4:1, resonance (lUHR; Schwarz 1984; 
Byrd et al. 1998). Athaneissoula et al. (1982) have discussed 
that for a flat rotation curve, where k — \/2fi, the ratio of 
the outer Lindblad, or 2:1, resonance (OLR) to lUHR radius 
is 

Rolr/Riuhr = ^"^^{^ = 2.64, (1) 
l-\/2/4 

while the ratio of the OLR to CR radius is 

Rolr/Rcr = 1 + V2/2 = 1.71. (2) 

This implies that the ratio of the CR to lUHR radius is 

Rcr/Riuhr = ^ r- , ^ 1-55. (3) 

l-^/2/4 

Since bars are encircled by inner rings, if present, Rbar must 
be similar to the inner ring radius. This suggests that TZ ~ 
1.6 or that bars are slow rotators in general. Statistics of 
observational data show that the median ratio of outer to 
inner ring diameter is 2.2 (Buta 1995) or the approximate 
mean of the theoretical OLR:IUHR and OLR:CR values. If 
outer rings in general can be linked to the OLR, this suggests 
a nearly even distribution of galaxies with fast or slow bars. 



Recent evidence suggests that outer rings may not always 
be linked to the OLR as expected though (Treuthardt et al. 
2008). 

Q,p, or TZ, can be estimated using a number of methods 
(e.g., Canzian 1993; Buta & Combes 1996; Puerari & Dottori 
1997; Salo et al. 1999; Weiner et al. 2001; Egusa et al. 2004; 
Zhang & Buta 2007) with the most direct being that of 
Tremaine & Weinberg (1984; TW). Much emphasis has been 
placed on this method recently but it has been suggested 
that morphological matching through dynamical simulations 
yield similar estimates while lacking the limited applicabilty 
of the TW method (Treuthardt et al. 2009). 

The goal of our study is to estimate the bar pattern 
speeds and measure the spiral arm pitch angles of 40 barred 
spiral galaxies. This will indicate which of these galaxies 
have low central dark halo densities and provide an estimate 
of the central SMBH mass, respectively. Based on the sim- 
ulation results of Booth & Schaye (2010), we expect to find 
that for a given halo mass, SMBHs in galaxies with fast bars 
will have a less-than-average mass due to anemic growth 
within the low central density dark halo. This means that 
the total luminous and dark mass of these galaxies (Mtot) 
predicted from the Mbh — Mtot relationship (Bandara et 
al. 2009) would be the minimum expected observationally 
measured value. We will investigate this claim with mea- 
surements of Mtot in a forthcoming paper. 



2 ANALYSIS 

2.1 Pitch Angle Measurements 

_B-band images from the Carnegie-Irvine Galaxy Survey 
(CGS; Ho et al. 2011) were used to determine P for our 
sample of barred spiral galaxies, since P has been shown to 
be independent of the wavelength at which it is measured 
(Seigar et al. 2006). These deep, high resolution images high- 
light the young stellar component of the spiral arms and of- 
fer a larger field of view than the i^s-band images, thereby 
revealing more of the outer spiral structure of each galaxy. 
The orientation and ellipticity values of the outer _B-band 
isophotes were determined with ELLIPSE in IRAF and used 
to derive an inclination (Hubble 1926; see Table 1). The im- 
ages were then deprojected to face-on by assuming that spi- 
ral galaxy disks are intrinsically circular. Additionally, the 
output from ELLIPSE was used to estimate Rbar in each 
galaxy. The radius, near the visible ends of the bar, where 
the ellipticity reaches a maximum was selected as Rbar- 

P was measured using two-dimensional fast Fourier de- 
compositions of the deprojected S-band images and assum- 
ing logarithmic spirals (Schroder et al. 1994; see Table 1). 
The Fourier fits were applied to visually selected annulus 
regions that range from just beyond the ends of the bars to 
the outer limits of the visible arms. P is then determined 
from peaks in the Fourier spectra as this is the most pow- 
erful method for finding periodicity in a distribution (Con- 
sidere & Athanassoula 1988; Garcia-Gomez & Athanassoula 
1993). The largest source of error in measuring P presum- 
ably comes from this selection of radial range (Seigar et al. 
2005, 2008), but galaxies with large inclinations (> 60°) also 
have variances in P of 10% (Block et al. 1999). For details 
on how accurate our method is for measuring P, see Davis 
et al. (2012). 
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Table 1. Estimated Galaxy Parameters. 

Explanation of columns: (1) Galaxy designation; (2) inclination in degrees derived from fitting ellipses to B- 
band isophotes; (3) bar semimajor axis length in the sky plane, estimated from ellipse fitting, and measured 
in arcseconds and converted to kpc (paranthetical values) using the default (i.e. Hq = 73.0, Omatter = 0.27, 
and ^vacuum = 0.73) cosmology-corrected scaling factors from NED; (4) D25 average isophotal diameters 
at 25 mag/arcsec^ in the B-band taken from HyperLeda (Paturel et al. 2003) and converted to kpc as in 
(3); (5) estimated value and la error of RcrI Rbar\ (6) spiral arm pitch angle measured by the authors 
of this paper, "Seigar et al. 2006, or by ''Block et al. 1999; (7) SMBH mass and mean error in 10^ solar 
masses derived from P using equation 2 of Seigar et al. 2008; (8) total mass of the galaxy and mean error 
in 10" solar masses derived from Mgu using equation 8 of Bandara et al. 2009 (*indicates an estimated 
minimum halo mass due to 7J. < 1.4 according to Booth & Schaye 2010); (9) de Vaucouleur morphological 
type index from HyperLeda (Paturel et al. 2003); (10) the maximum of Qx in the bar region estimated from 
mass models derived from its-band images. 



2.2 Determining the Potential 

Attempting to recreate the observed morphologies of our 
sample of galaxies through simulation modeling requires us 
to first derive gravitiatonal potentials. Assuming that the 
near-IR light distribution follows the mass distribution, we 
applied a mass-to-light ratio [M/L) correction to the CGS 
near-IR Ks-haxiA images of our galaxies to obtain surface 



mass densities. The M/L correction was estimated from 
the azimuthally averaged CGS B — V color index of each 
galaxy and applying the K-haxiA correction given by Bell et 
al. (2003). 

Next, a two-dimensional iterative bulge-bar-disk decom- 
position program (see Laurikainen et al. 2004 for a descrip- 
tion) was applied in order to more accurately model the 
bulge component of the galaxies. The bulge component was 



© 2002 RAS, MNRAS OOO.mfTTi 



4 P. Treuthardt et al. 



then removed, each galaxy was deprojected to a face-on ori- 
entation, and the bulge was then added back in. We assume 
the bulges to be spherical, so this process eliminates any 
bulge stretching from the deprojection process. 

The deprojected image was then approximated by the 
Fourier decomposition method described by Laurikainen & 
Salo (2002) and the disk gravity was derived from the even 
m = to 8 components. The empirically derived ratio of 
vertical to radial scalelength corresponding to morphological 
type (de Grijs 1998) was used to apply a constant scaleheight 
throughout the disk. Our mass models also include a dark 
halo component based on the universal rotation curve of 
Persic et al. (1996). The halo profile was calculated using 
L/L* derived from HyperLeda (Paturel et al. 2003) data and 
the average distances calculated from NED (NASA/IPAC 
Extragalactic Database) data. 

A limited amount of kinematic data exists for our sam- 
ple of galaxies (Mathewson et al. 1992; Marquez et al. 2004; 
Catinella et al. 2005). What data we have are chiefly in 
the form of Ha rotation curves which only extend to about 
80% of the visible galaxy disk (Mathewson et al. 1992) and 
do not extend into the dark matter dominated regions of 
our mass models. The available data, therefore, does little 
to constrain the outer mass-modeled circular speed curves. 
The upper plot of Figure 1 shows an example of azimuthally 
averaged circular speed curve derived from the mass model 
of ESO 506-004. The mass model was scaled so that the 
circular speed curve would more accurately fit the rotation 
curve data derived from Figure 3 of Mathewson et al. (1992). 
The lower plot in Figure 1 shows the Lindblad precession 
frequency curves derived from the circular speed curve. Our 
estimate of Q.p, average error, and bar semimajor axis radius 
for this galax;y are overlayed. The radial regions where spe- 
cific resonances in the linear (epicyclic) approximation are 
located are indicated where the horizontal line correspond- 
ing to Q.p intersects each frequency curve. For this galaxy, 
we estimate that the bar extends to the lUHR and not to 
the CR. 

Once gravitational potentials were determined for our 
sample, we were also able to estimate the strength of the bar 
in each galaxy. This was done by calculating two dimensional 
maps of the radial (Fr) and tangential {Ft) forces. The 
tangential forces induced by the bar were calculated using a 
Polar method discussed by Laurikainen & Salo (2002) and 
Laurikainen et al. (2004). The Polar method was applied to 
an azimuthal Fourier decomposition of the intensity and the 
even components, up to m = 20, were then converted to the 
corresponding potential components (Salo et al. 1999). The 
ratio of the maximum of Ft to the azimuthally averaged 
Fh at each radius produces a radial profile of the distance 
dependent maximum tangential force. The maximum of the 
this profile gives a single measure of the bar strength, Qg. 
Qg, however, can be affected by spiral arm torques, making 
it a convolution of the bar (Qb) and spiral strength (Qs; 
Buta et al. 2005). Therefore, Qg can be thought of, at worst, 
as an upper limit to the true bar strength, Qb- 

2.3 Simulations 

Our sample of galaxies was modeled by simulating the be- 
haviour of a two-dimensional disk of inelastically colliding 
(or "sticky"), non-self-gravitating gas particles in our pre- 



determined potentials. The details of the simulation code we 
used can be found in Salo et al. (1999) and Salo (1991). The 
advantage of using a sticky particle code over a smoothed 
particle hydrodynamic (SPH) code, for example, is that it 
better approximates the cold gas component of a galaxy 
(Merlin & Chiosi 2007). In SPH, gas particles are subject to 
nongravitational forces, such as pressure gradients, making 
it better suited for modeling the warm gas phase. Therefore, 
it is appropriate to use sticky particle models for matching 
observed morphological components containing star forma- 
tion, such as resonance rings and spiral arms. The overall 
modeling process closely followed that of Treuthardt et al. 
(2009). For example, we produced snapshots of 100,000 gas 
particles at the bar rotation period of interest by simulating 
the behaviour of a two-dimensional disk of 20,000 gas parti- 
cles and aligning and coadding individual snapshots within 
± 0.2 bar rotation periods in 0.1 period increments. We also 
assumed that each galaxy has only one pattern speed, that 
of the bar, and Qp was the main parameter that was varied. 
With an estimate of Rtar and using the Lindblad precession 
frequency curves derived from our potentials, we are able to 
link Hp to TZ via Rcr- In doing so, we varied Qp by allow- 
ing TZ to increase from 1.0 by increments of 0.1. Treuthardt 
et al. (2008) investigated the morphological effects of dif- 
fering the gravitational potential used in the simulations of 
NGC 1433. They found that the different M/L models of 
Bell & de Jong (2001) had little effect on the modeled mor- 
phology of the galaxy. The amplitude of the bar potential 
was also incrementally altered in order to account for possi- 
ble uncertainties in the bar height or halo contribution. The 
morphology of the resultant models were markedly differ- 
ent from each other in terms of inner ring shape and outer 
spiral morphology (see Figure 9 of Treuthardt et al. 2008). 
This implies that the use of grossly inaccurate gravitational 
potentials will not produce models with morphologies that 
can be easily matched to observations. Further examples of 
morphological variation within sticky particle models of a 
galaxy can be found in Treuthardt et al. (2008). 

The best fitting TZ of each galaxy was determined by vi- 
sually comparing the observed B-band morphology to that 
produced in a series of models. Deep B-band images with 
good resolution are ideal for comparing observed low ve- 
locity dispersion, young stellar components, such as reso- 
nance rings, to our models. In our simulations, the viscosity 
within various regions of the gas particle system may not 
be comparable to real systems, meaning features seen in 
the simulation snapshots can correspond to different phys- 
ical timescales in different regions. With that in mind, we 
primarily compared the morphology in the region near the 
ends of the bar (such as the size and shape of any inner 
rings) but also considered the outer spiral structure, espe- 
cially any outer ring features. Since we are dealing with gas 
particle simulations of barred spiral galaxies, areas within 
the bar region of our models tend to be void of gas parti- 
cles and are not considered in matching with the observed 
images. Four independent estimates of the model morphol- 
ogy that best matches the observed morphology were used 
to determine the weighted average value and la error of TZ 
given in Table 1. In some cases a slower Qp provided a better 
fit to the outer spiral structure while a faster Qp provided 
a better fit to the inner structure. This is not unexpected 
since N-body simulations have shown that the spiral com- 
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Figure 1. Plots of the rotation curve (top) and Lindblad precession frequency curves (bottom) for ESO 506-004. Using the same scaling 
as in column 3 of Table 1, 10 arcseconds corresponds to 2.85 kpc. The upper plot shows the total, bulge, disk, and dark matter halo 
components of the azimuthally averaged mass modeled circular speed. The dark matter halo component was derived from the universal 
rotation curve model of Persic et al. (1996). Ha rotation curve points derived from Fig. 3 of Mathewson et al. (1992) are overlayed with 
the estimated errors in their values. The inner region of this specific galaxy appears to have little hydrogen gas. The modeled major 
axis velocity profile (points), derived from the modeled gas particle distribution of our best simulation, is also overlayed. The velocity 
profile data were extracted from within 3° of the major axis position angle and corrected for inclination. The modeled velocity profile 
appears to agree with the observed rotation curve. The lower plot shows Q (the circular frequency), Q ± k/2 (where k is the epicyclic 
frequency), Q ± k/4, and the bar pattern speed {Qp = 44.0 km s"'^ kpc"'^) corresponding to our estimate of TZ = 1.74. The position of 
the short vertical line associated with Qp indicates the estimated bar semimajor axis length, while the length of the line corresponds to 
the average error in Qp. For this galaxy, it appears the bar semimajor axis extends to the inner 4:1 resonance. 



ponent may have a slower pattern speed then that of the 
bar (e.g. Rautiainen & Salo 1999). This may be the reason 
why we tend to find few galaxies in our sample with con- 
vincingly fast bars. Figure 2 shows the deprojected B-band 
and best model images corresponding to our estimated TZ. 
Subsection 2.3.1 provides brief discussions on the individual 
galaxies and the corresponding models. 

A number of our sample galaxies appear to have dust 
lanes along the bar major axis when examined in the 
high resolution B-band images from CGS. We make a 
note of these dust lanes below in the comments on in- 
dividual galaxies. Athanassoula (1992) has experimented 
with two-dimensional time- dependent hydrodynamical sim- 
ulations that produce dust lanes in barred spiral galaxies. 
The gravitational potentials used in the simulations consist 
of two axisymmetric components, which approximate the 
bulge and disk potential, and a Ferrers ellipsoid, which ap- 
proximates the bar potential. A dark halo component was 
not included in the models. Based on the dust lane mor- 
phology observed in her models, Athanassoula concludes 
that bars containing dust lanes must be in the fast domain 
{TZ = 1.2 ± 0.2). She claims the convex curvature of dust 
lanes shown in her slow bar models are not observed in real 
galaxies and concludes that slow bars do not exist. The dif- 
fuse nature of the few dust lanes seen in our sample of galax- 
ies generally makes the morphology of the dust lanes difficult 
to determine visually. 



2.3.1 Comments on Individual Galaxies 

ESO 121-026: Our best models were not able to repro- 
duce the outer flocculent arms of this SB(rs)bc galaxy (de 
Vaucouleurs et al. 1991). The inner ring, however, is well 
matched by the model shown in Figure 2. The model also 
displays a pointy ring elongated perpendicularly to the ma- 
jor axis and encircling the inner ring. This pointy ring is also 
seen in the observed image. The B-band image displays a 
possible diffuse dust lane extending along only the bottom 
half of the bar's major axis. 

ESO 380-001 : The inner pseudoring and faint outer ring 
of this galaxy, seen in the B-band image, is reflected in the 
(R')SB(s)b classification by de Vaucouleurs et al. (1991). 
In the observed image, the spiral arms can be traced from 
their origin at the ends of the bar, through their clockwise 
winding around the bar, and finally transitioning into the 
outer ring. Our models display a similar morphology but 
also include an obvious inner ring. This may be hinted at in 
the B-band near the bar end closest to the bottom of the 
image. Symmetric dust lanes are also seen extending along 
the major axis of the bar in the B-band image. 

ESO 506-004: The model shown for this SAB(r)ab 
galaxy (de Vaucouleurs et al. 1991) displays a good match to 
the observed inner ring and the strong spiral arm originating 
at the left end of the inner ring major axis. 

IC 1953: Our models of this SB(rs)d galaxy (de Vau- 
couleurs et al. 1991) produce a series of nested rings sur- 
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Figure 2. Deprojected _B-band and best gas particle model images of our sample of galaxies. The models shown correspond to our 
estimated values of TZ.. Some B-band images display black masks that were used when a foreground star could not be otherwise removed. 
The thick tick-marks along the right edge of the gas particle images indicate the diameter of corotation. 



rounded by outer flocculent arms. The outermost modeled 
ring is a hexagonal structure whose major axis appears ro- 
tated slightly counter-clockwise to the bar major axis. Hints 
of this outer most ring structure along the major axis is seen 
in the B-band image. Within this elongated hexagonal ring 
is a nearly circular ring surrounding a ring elongated along 
the bar major axis. Hints of this nearly circular ring, but not 
the inner elongated ring, are also present in the observed im- 
age. 

IC 2367: The size and extent of the inner ring of this 
SB(r)b galaxy (de Vaucouleurs et al. 1991) is well matched 
by the model shown. 



IC 2560: The inner ring of this (R')SB(r)b galaxy (de 
Vaucouleurs et al. 1991) is well reproduced by our best es- 
timate model. Also, there appear to be dust lanes running 
along the major axis of the bar when viewed in the B-band. 

IC 5240 : Buta et al. (2007) classify this galaxy as 
SB(r)a with a clear inner ring and flocculent exterior spi- 
ral structure. Upon closer inspection, the upper half of the 
inner ring appears to consist of three straight arm segments. 
Our simulations fail to produce models that directly recreate 
the strong, nearly circular inner ring seen in the deprojected 
B-band image. Our best model does approximate the in- 
ner ring morphology from a series of disjointed spiral arms. 
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There appears to be a dust lane along only the right half 
of the bar's major axis which extends into dust in the inner 
ring. 

IC 5273: The B-band image of this SB(rs)cd galaxy 
(de Vaucouleurs et al. 1991) does not display any obvious 
resonance features such as rings. We attempted to match 
the inner spiral arms, near the bar semimajor axis length, 
of our models to that in the observed image. For example, 
the observed spiral arm originating near the bottom end of 
the bar major axis is matched in our given model. 

NGC 0151: Our models are not able to recreate the 
asymmetry observed in the outer spiral structure of this 
SB(rs)bc galaxy (Buta et al. 2007) since only the even 
Fourier m values were used to calculate the disk gravity. 
This causes our models to be bisymmetric. Our estimate of 



TZ came from matching the the observed and modeled inner 
ring morphology. 

NGC 0337: This SB(s)d galaxy (de Vaucouleurs et al. 
1991) has the strongest bar of our sample {Qg = 0.800). Our 
models produce a small inner ring not visible in the B-band 
image. Our estimates of TZ were based on comparing the 
modeled 1/4 and L5 Lagrangian regions and the morphology 
of the spiral arms extending from the inner ring to that in 
the observed image. 

NGC 0782: This SB(r)b galaxy (de Vaucouleurs et al. 
1991) appears to be a three armed system when observed 
in the _B-band. As was discussed regarding NGC 0151, our 
models are necessarily bisymmetric and will not produce 
three armed systems. Our best estimate of TZ came from 
attempting to match the inner rings and the arm structure 
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originating at the bottom of the inner ring in both the ob- 
served and modeled image. 

NGC 0945: de Vaucouleurs et al. (1991) classify this 
galaxy as SB(rs)c. Our models are able to reproduce the 
shape and size of the observed inner ring but the outer spiral 
structure is more open than seen in the observed image. This 
may be the result of NGC 0948 which is visually nearby 
and classified as having an uncertain SB(s)c morphology (de 
Vaucouleurs et al. 1991). NGC 0945 also appears to have 
asymmetrical dust lanes running along the length of the bar 
major axis. 

NGC 1022: Buta et al. (2007) classify this galaxy as a 
peculiar (R)SB(rs)a. Our model, corresponding to our best 
estimate of TZ, provides a good match to both the inner and 
outer ring morphology. 

NGC 1317: Buta et al. (2007) describe this galaxy 
as having a double-barred morphology and classify it as 
(R')SAB(rl)a. We attempted to recreate the observed mor- 
phology by varying Q.p of the nuclear bar, which is the weak- 
est in our sample (Qg = 0.040), rather than the less promi- 
nent broad oval, or primary bar. Our models produced a 
tight spiral structure around the nuclear bar with a higher 
density of gas particles on the leading edge of the spiral 
and a lower density on the trailing edge. We matched the 
lower density trailing edge of our models to the gaps seen 
in the B-band emission surrounding the nuclear bar region. 
Sticky-particle modeling of this galaxy by Rautiainen et al. 
(2008) led to an estimate of 7^ = 0.89 ± 0.25, nearly half of 
our estimate. Although Rautiainen et al. (2008) do not com- 
ment on this specific galaxy, it is apparent from Figure 4 of 
their paper that their models were produced by varying Q,p 
of the broad oval rather than the nuclear bar. The nuclear 
bar does not display any apparent dust lanes, but a nearly 
symmetrical dust spiral is seen in the broad oval region. 

NGC 1723: This peculiar SB(r)a galaxy (de Vau- 
couleurs et al. 1991) may be the product of an interaction 
with a nearby galaxy. The fact that one of the outer spiral 
arms extends and curves toward a smaller field galaxy in our 
raw B-band image is apparent evidence for this. Our sim- 
ulations were not able to recreate the observed outer spiral 
structure or the slightly oblong nature of the inner ring seen 
in our deprojected image. The model corresponding to our 
weighted mean value of TZ contains an inner ring of similar 
size to the observed inner ring. Larger or smaller values of 
TZ result in models with larger or smaller inner rings, respec- 
tively. 

NGC 1832: de Vaucouleurs et al. (1991) classify 
this galaxy as SB(r)bc. Our model, corresponding to our 
weighted mean estimate of TZ, produces an inner ring of 
similar extent to that seen in the observed image. The inner 
ring in the B-band image may be the result of gas circulat- 
ing along so called "banana" orbits around the Z/4 and 1/5 
Lagrangian points. This is also seen in our model, but the 
density of gas particles is not very large in the portions of 
these orbits that bring them close to the bar. Rautiainen et 
al. (2008) estimate TZ to be 1.74 ± 0.40 from their sticky- 
particle simulations using the same code as the current sim- 
ulations. Within the errors, this coincides with our estimate 
of7^ = 1.61 ±0.08. 

NGC 2223: Our model of this SAB(r)b galaxy (de Vau- 
couleurs et al. 1991) produces a strong inner ring of the same 
size and extent as the one that is hinted at in B-band im- 



age. Our model also recreates the inner spiral arm structure 
observed just beyond the ends of the inner ring major axis. 

NGC 2525: This SB(s)c galaxy (de Vaucouleurs et al. 
1991) is well matched by our model. The size and position 
of the modeled LA and L5 Lagrange regions duplicate that 
seen in the B-band image. The modeled spiral structure also 
closely matches the observations, although the gas particle 
density is not high. 

NGC 2763: de Vaucouleurs et al. (1991) classify this 
galaxy as a peculiar SB(r)cd. The simulation model of this 
galaxy, shown in Figure 2, is the weighted mean of individ- 
ually estimated TZ values ranging from 1.2 ±0.2 to 1.9 ±0.2. 

NGC 3124: Buta et al. (2007) classify this galaxy as 
SB(r)bc. This galaxy does not host a traditional bar but 
rather a bar-spiral hybrid winding in the opposite sense as 
the outer spiral arms. This is especially apparent in the Ks- 
band image. Our models do not reproduce the nearly circular 
looking inner ring seen in the B-band but rather approxi- 
mate it as a pair of disjointed spiral arms. A one-armed dust 
spiral is seen in the B-band winding in the same sense as 
the stellar spiral arms. This one-armed dust spiral is seen 
crossing the left semi-major axis of the bar. 

NGC 3275: The apparent inner ring seen in the B- 
band image of this SB(r)ab galaxy (Buta et al. 2007) is 
recreated by our best estimate model in the form of spiral 
arm segments outside of a true inner ring. A faint hint of 
a thin "true" inner ring is seen in the observed image near 
the left end of the bar, inside the apparent inner ring. A 
morphological match is also seen when comparing the outer 
spiral arm that originates from near the right end of the 
bar. Models produced by Rautiainen et al. (2008) estimate 
TZ — 1.53 ± 0.34 which coincides with our estimated value 
within the errors. The B-band image of this galaxy shows 
diffuse, nearly symmetrical dust lanes and spirals in the bar 
region. The dust lanes are originate near the ends of the bar 
and transition into spirals at the bar ends. 

NGC 3347: Our assumption of intrinsically circular 
galaxy disks does not appear to be valid for this galaxy 
which leads to a somewhat unrealistic deprojected image. 
NGC 3347 is a member of a group (Garcia 1993) and ap- 
pears, in the CCS images, to be interacting with a nearby 
peculiar spiral galaxy of uncertain morphology (NGC 3354; 
de Vaucouleurs et al. 1991). Our simulation models were not 
able to recreate the tightly wrapped outer spiral arms of this 
SB(rs)b galaxy (de Vaucouleurs et al. 1991). The shape and 
extent of the modeled inner ring closely matches that sur- 
rounding the observed bar. 

NGC 3450: The observed morphology of this SB(rs)b 
galaxy (Buta et al. 2007) is well matched by our model. The 
inner ring, as well as the outer spiral structure, appear to 
agree in both the observed and modeled image. A dust lane 
elongated nearly perpendicular to the semi-major axis of the 
lower left end of the bar is seen in the B-band image. 

NGC 3513: This galaxy is classified as SB(s)c by Buta 
et al. (2007) and contains one of the strongest bars of our 
sample {Qg = 0.540; see Table 1). Our simulations were un- 
able to easily reproduce the main spiral arms alone. Instead, 
most of our gas particle models show diffuse arms emanat- 
ing perpendicularly to the bar from an apparent nuclear 
ring. In the model associated with our weighted average es- 
timate of TZ, a smaller number of gas particles also appear 
to emanate from this ring, run parallel to the bar, curve 
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around the bar ends, and finally reconnect to the stronger, 
more diffuse arms on the opposite side of the ring. This less 
prominant gas particle distribution more closely resembles 
the observed morphology than the more obvious gas parti- 
cle distribution. Sticky-particle modeling by Rautiainen et 
al. (2008) produces an estimate of 7?. = 1.50±0.29. Our esti- 
mate of 7?. = 1.24 ± 0.07 falls within the error of Rautiainen 
et al.'s estimate. 

NGC 3660: Recreating the outer spiral arm morphol- 
ogy of this SB(r)bc galaxy (Buta et al. 2007) proved to be 
difficult. Our weighted mean estimate of fip, or TZ, provides 
a good morphological match of the inner ring between both 
the observation and the model. The B-band image shows 
faint, symmetrical, and diffuse dust lanes running along the 
major axis of the bar. 

NGC 3881: Our models of this SB(rs)bc galaxy (Buta 
et al. 2007) produce an inner ring that is not visible in the 
B-band image. The extent the inner ring shown in our model 
does correspond to the extent of the straight portion of the 
symmetrical dust lanes seen along the major axis of the bar. 
The turn-over radius of these dust lanes is similar to the 
inner ring semimajor axis length. The outer spiral structure 
also approximates the observed spiral structure near the bar. 

NGC 4050: The morphology of this SB(r)ab galaxy (de 
Vaucouleurs et al. 1991) is well matched by our model. Both 
the inner ring and the outer spiral structure appear to agree. 
A single diffuse dust lane is seen in the B-band running from 
the nucleus along the right bar semi-major axis. 

NGC 4593: Buta et al. (2007) classify this galaxy as 
(R')SB(rs)ab. Our model does well to approximate the shape 
and extent of the inner and outer ring. The outer ring ap- 
pears to be due to gas trapped in "banana" orbits and cir- 
culating around the 1/4 and L5 Lagrangian regions, as was 
seen in the simulation models of NGC 1433 (Treuthardt et 
al. 2008). The observed image of NGC 4593 also shows an 
apparent "plume", similar to those seen in NGC 1433, near 
the leading edge of the bar and in the upper left quadrant of 
the _B-band image. Our model also produces a linear compo- 
nent to the "banana" orbits, close to the outer boundary of 
the inner ring. This feature is also seen in lower half of the 
observed image. A strong dust lane is seen along the upper 
left semi-major axis of the bar, while a fainter, more diffuse 
dust lane is seen along the opposite semi-major axis. 

NGC 5135: This SB(s)ab galaxy (de Vaucouleurs et al. 
1991) was well matched by our model. Both the inner ring 
and outer spiral morphology are reproduced. The inner re- 
gions of this galaxy display a lot of dust in the _B-band, but 
no dust lanes are apparent. 

NGC 5156: We were not able to reproduce the over- 
all outer spiral structure observed in this SB(r)b galaxy (de 
Vaucouleurs et al. 1991). Determining a good morphological 
match to the inner ring also proved to be difficult as our 
estimates oiTZ ranged from 1.4±0.2 to 1.8±0.2. The model 
we show, corresponding to the weighted mean of our TZ val- 
ues, does approximate the two observed spiral arms closest 
to the bar, in the lower half of the image. 

NGC 5339: This peculiar SB(rs)a galaxy (de Vau- 
couleurs et al. 1991) is well matched by our model. The 
outer spiral arm morphology is recreated by a series of arms 
that form a slight hexagonal pattern. The extent and shape 
of the inner ring is also approximated. Nearly symmetrical 



dust lanes are seen in the B-band along the major axis of 
the bar. 

NGC 5728: The inner ring of this (Ri)SAB(r)a galaxy 
(Buta et al. 2007) is reproduced by our model. An outer 
ring is also created by our simulation, but the observed ring 
is much more circular in the deprojected image. Symmetri- 
cal dust lanes are also seen along the bar major axis when 
viewed in the B-band. 

NGC 5938: A large number of foreground stars are in 
the field of this SB(rs)bc galaxy (de Vaucouleurs et al. 1991). 
Those stars that could not be blended to the background 
where masked, leaving obvious artifacts in the image. The 
spiral structure far from the bar is not well matched by our 
model. The left side of the observed apparent inner ring 
is matched by our modeled inner ring. The right side of 
the observed apparent inner ring is better matched by the 
modeled spiral arm originating just outside the inner ring. 

NGC 6782: This galaxy is classified as (RiR2)SB(r)a 
by Buta et al. (2007). The outer ring features are best fit by 
7?. ~ 1.3 while the inner ring is best fit by 7?. ~ 1.1. In either 
case, the bar is still classified as fast. For this galaxy, Rauti- 
ainen et al. (2008) estimate TZ = 1.25 ± 0.23 but they were 
not able to recreate the outer pseudo-ring in their models. 
Our models display this outer pseudo-ring with a very sim- 
ilar value of TZ. Lin et al. (2008) produced two-dimensional 
hydrodynamical simulations of this galaxy and were able to 
recreate the symmetric dust lanes observed in the B-band. 
They, however, found 7?. to be 1.5. 

NGC 6923: Our model morphology matches the ob- 
served inner ring and lower spiral arm of this SB(rs)b galaxy 
(de Vaucouleurs et al. 1991). 

NGC 7059: This highly inclined {i = 71.1°; see Table 1) 
SAB(rs)c (de Vaucouleurs et al. 1991) galaxy displays little 
structure in its deprojected B-band morphology that could 
be thought of as an unambiguous, organized product of flp. 
All of our simulations produce models with a clear inner ring 
which is not evident in the B-band. Many of our models also 
display four low particle density spiral arms extending from 
near the ends of the major axis of the inner ring. This spiral 
structure may, debatably, be hinted at in the deprojected B- 
band image. Consequently, the weighted mean of our best 
morphological matching estimates {TZ ~ 1.7) was deter- 
mined from individual estimates ranging from TZ — 1.9 ±0.1 
to 7?.= 1.2 ±0.2. 

NGC 7070: de Vaucouleurs et al. (1991) classifies this 
galaxy as a non-barred SA(s)cd. Both the B and Ks- 
band images show a clear, highly elliptical structure at the 
galaxy's center. We attempted to match the observed spi- 
ral arm morphology, but this galaxy's flocculent structure 
and lack of any obvious resonance features induced a large 
amount of error in our TZ estimates. The errors of our indi- 
vidual estimates ranged from ~ 12% — 33%. 

NGC 7218: The morphology of the pair of spiral arms 
originating from the left end of the observed bar of this 
SB(rs)cd galaxy (de Vaucouleurs et al. 1991) is approxi- 
mated by our model. The more tightly wound of these arms 
appears to connect to a Lagrange region, and continue on 
to form a wider, more diffuse, and more open spiral arm in 
our model. The origin point of this wider arm may be seen 
near the upper right end of the bar in the B-band image. 

NGC 7392: A central oval component of this SA(s)bc 
galaxy (de Vaucouleurs et al. 1991) is more apparent in the 
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Ks-hand image than in the _B-band. Our model does well 
to match the inner pseudo ring, as well as the outer spiral 
structure. Particularly, our model provides a morphological 
match to a looping ring seen in the upper half of the galaxy 
image. Faint, asymmetrical dust lanes are seen in the B- 
band along the major axis of the central oval. 

NGC 7723: Buta et al. (2007) classify this galaxy as 
SB(rs)b. Our model is able to recreate the inner ring and 
outer spiral structure from a series of complete and seg- 
mented spiral arms. Modeling by Rautiainen et al. (2008) 
produce an estimate of TZ = 1.31±0.29 which coincides 
with our model-based estimate within the errors. Dust lanes 
are seen in the _B-band image of this galaxy. The dust lane 
that runs along the leading edge of the right semi-major 
bar axis is quite apparent and linear. The dust lane running 
along the opposite end of the bar appears to originate from 
the trailing side of the bar and cross to the leading side. 



3 RESULTS AND DISCUSSION 

Booth & Schaye (2010) have developed models of black hole 
growth in galaxies that appear to reproduce the empirical 
relationship found between black hole mass and halo mass 
(Bandara et al. 2009; see figure 3 of Booth & Schaye 2010). 
The implication from their models is that galaxies with low 
dark halo concentrations will host SMBHs with a range of 
masses, but they will be less massive than average for a given 
halo mass (see Figure 3). 

A simple, easily applicable tool for determining the 
SMBH mass of nearby spiral galaxies is through the rela- 
tionship between Mbh and P found by Seigar et al. (2008). 
The Mbh — P relationship is fit by Seigar et al. (2008) using 
a double-power-law model given as 

Mbh = 2^'-^^'- Mbh, (^) [l + (^) ] , (4) 

where /3 is the slope of the power law for large pitch angles, 
7 is the slope of the power law for small pitch angles, Pb 
is the transition from small to large pitch angles, a governs 
the sharpness of the transition, and Mbh, is the black hole 
mass for a pitch angle Pj. The values of the parameters that 
provide this best fit model are /3 = 126.1, 7 = 2.92, Ph = 
40.8°, a = 23.5, and Mbh, = 1.72 x lO^M© (see equation 
2 of Seigar et al. 2008). Other model fits have also been 
suggested such as a linear (Seigar et al. 2008) or cotangent 
relation (Ringermacher & Mead 2009). 

In Figure 4 (left) we plot our estimated values of TZ ver- 
sus our measurements of pitch angle. There does not appear 
to be any correlation between either fast, intermediate, or 
slow bar pattern speed and pitch angle. Those galaxies with 
clearly fast bars span a range in P from 13.4° for NGC 6782 
to 32.9° for NGC 0337. The two galaxies with the smallest 
pitch angles of those with fast bars are NGC 5135 and NGC 
6782. The outer spiral arms of these galaxies are not log- 
arithmic as we assume and appear to be the result of gas 
trapped in resonance regions associated with Sip (see Figure 
5). The relative amplitudes of the m = 2 and m = 4 Fourier 
components indicate that a single spiral pattern is dominate 
(see Figure 6) meaning that P is not constant with radius in 
our measurements of these two galaxies. When pitch angle 
is translated to SMBH mass via equation 2 in Seigar et al. 



(2008) we can see that Mbh also does not correlate with 
bar pattern speed (Figure 4, right). That is, galaxies with 
fast bars, and low central density dark matter halos (Debat- 
tista & Sellwood 2000), have SMBH masses similar to other 
barred spiral galaxies. 

If the assertion is true that galaxies with low dark halo 
concentrations host black holes less massive than average for 
a given halo mass and the Mbh — P relationship holds, then 
we can make a prediction of the minimum mass of our fast- 
bar galaxies. Figure 7 shows our TZ estimates versus the total 
galaxy mass implied from the Mbh — P (Seigar et al. 2008) 
and Mbh — Mtot (Bandara et al. 2009) relationships. Here 
Mtot is the total gravitational mass of the galaxy within 
the radius where the density profile of the galaxy exceeds 
the critical density of the universe by a factor of 200. The 
predicted total mass of those galaxies with convincingly fast 
bars is expected to be a minimum (see Table 1). Follow-up 
rotation curve measurements will test this prediction. 

A plot of P versus the de Vaucouleurs morphological 
type index (T; see Figure 8) shows a trend in agreement 
with Kennicutt (1981) and Seigar & James (1998b). Gala:x- 
ies identified as early-type have correspondingly small values 
of P while late-types have larger values of P. The trend be- 
comes a tighter correlation when the focus is shifted from the 
overall sample to only those galaxies with clearly fast bars. 
A linear fit to the data of Roberts et al. (1975) concerning 
theoretical pitch angle versus Hubble type is also overlayed 
(dot-dashed line). The y- intercept of this fit was systemati- 
cally lower than our data points, so it was multiplied by an 
ad hoc factor of ~ -0.46. The slope of this theoretical fit is 
more similar to the slope of the linear fit of the galaxies with 
fast bars (dashed green line) than to the slope of the linear 
fit of the entire sample (dotted line) . 

Das et al. (2008) have shown that a trend exists where 
bar strength (Qg) is anticorrelated with the central velocity 
dispersion (uc) of the bulge. This result suggests that the 
growth of a central mass may be closely linked to the evolu- 
tion of the bar. Additionally, Seigar et al. (2008) show that 
a good anticorrelation exists between CTc and P. This is a 
result of CTc also being an indicator of Mbh (Gebhardt et al. 
2000; Ferrarese & Merritt 2000). These results imply that a 
correlation between bar strength and P should be expected. 
In Figure 9 we plot Qg as a function of P and do not see 
this correlation in the overall sample. However, when focus- 
ing only on the small sample of galaxies with clearly fast 
bars, Qg appears to increase with increasing P. 

An overall trend becomes more apparent when we plot 
Qb as a function of T for a larger sample of galaxies (see Fig- 
ure 10). We have already shown that there is a loose correla- 
tion between P and T in general and perhaps a tighter cor- 
relation when only considering galaxies with fast bars. Also, 
if Qs is small, then Qg approaches Qb, effectively making 
Qg an upper limit of Qs. The larger overall sample displays 
a very weak correlation which agrees with the results of Lau- 
rikainen et al. (2007) and Seigar & James (1998a), although 
the latter authors determined the bar strength through an 
alternate method and plotted it as a function of bulge-to- 
disk ratio. When focusing once again on galaxies with clearly 
fast bars, it becomes more evident that a trend between bar 
strength and T, or P, exists for those galaxies. This seems 
to suggest that galaxies with a smaller dark matter halo 
contribution in the central regions follow a more predictable 
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Figure 3. Sketches of Figure 3 from Booth & Schaye (2010). The red hne represents the empirical SMBH and halo mass relationship 
(Bandara et al. 2009) indicating the average SMBH mass for a given halo mass. The green ellipse outlines the region of less-than-average 
SMBH masses due to less concentrated haloes. The upper figure shows that for a galaxy with a low concentration halo and a given 
measured halo mass, the SMBH mass is overestimated (A) when only the empirical relationship is used. In reality, a galaxy with a low 
concentration halo will have a lower SMBH mass (B) for the same measured halo mass. The lower figure shows that when the SMBH 
mass is measured for a galaxy with a low concentration halo, the empirical relationship implies a halo mass that is lower (A) than the 
actual halo mass (B). 



RcR/RbarVS. Pitch Angle Rcp/RtarVS. Supermassive BH Mass 
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Figure 4. Plots of estimated TZ = RcrI Rbar versus measured spiral arm pitch angle (left) and SMBH mass (right) for galaxies with 
fast, intermediate, or slow bar pattern speeds. The red open circles correspond to galaxies where TZ is estimated to be in the slow bar 
domain within the errors, the green filled circles correspond to galaxies in the fast bar domain within the errors, and the blue filled 
squares correspond to galaxies that could fall in either domain given the errors. The fit given by equation 2 in Seigar et al. (2008) was 
applied to the pitch angle measurements in the left plot to determine Mbh in the right plot. The dashed line indicates 7^ = 1.4, the 
boundry between fast and slow bars. There appears to be no correlation between bar pattern speed and pitch angle or black hole mass. 
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Figure 5. _B-band images of NGC 5135 (left) and NGC 6782 (right) with contours of the m=2 Fourier components overlayed. The 
contours are fit to the outer spiral arms with the assumption that they are logarithmic spirals. The figures show that the outer arms 
appear logarithmic at small radii but deviate from this assumption at larger radii. This is likely due to gas being trapped at the OLR 
and implies that P is not constant with radius in our measurements. 




trend in bar strength and spiral arm morphology than those 
with more dominant dark matter contributions. 

Finally, Figure 11 displays TZ versus T for our data as 
well as that of Rautiainen et al. (2008). The general trend 
of TZ increasing with an increasing T, seen in Figure 6 of 



Rautiainen et al. (2008) , is not as clear with the addition of 
our data. Instead, the few TZ values greater than ~2.0, seen 
from 3 < T < 5, appear to be anomalous rather than part of 
a general trend. The large range in T, or P, seen for slow, 
intermediate, and fast bar galaxies in this larger overall sam- 
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Figure 7. Plot of TZ = Rcr/ Rbar versus total galaxy mass showing the distribution of total masses for galaxies. The horizontal dashed 
line corresponds to 7?. = 1.4, or the boundary between the fast and slow bar domains 
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Pitch Angle vs. de Vaucouleurs Type Index 
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Figure 8. A plot of the measured spiral arm pitch angle versus de Vaucouleurs morphological type index taken from HyperLeda. The 
dashed green line is a fit through the points corresponding to galaxies with fast bars, the dotted line is a fit through all the points, and the 
dot-dashed line is the fit predicted by Roberts et al. (1975). The former two fits were obtained with the Interactive Data Language (IDL) 
program FITEXY. The latter, theoretical fit was systematically lower than the other two observationally derived fits. The y-intercept of 
the dot-dashed line was, therefore, scaled by a factor of -0.46 so that all three lines intersect at T 2.9. 



pie is indicative of the scatter around the empirical SMBH 
and halo mass relation of Bandara et al. (2009) predicted by 
Booth & Schaye (2010). 



4 CONCLUSIONS 

We set out to estimate the SMBH mass of a sample of 40 
barred spiral galaxies and determine which of these galaxies 
have low central dark halo densities. We did this by using 
a 2D FFT to measure P and estimated Mbh via the best 
model fit of the Mbh — P relationship (Seigar et al. 2008). 
We also created simulation models of these galaxies with 
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Figure 9. A plot of bar strength, Qg, versus pitch angle. The anticorrelation of Qg versus CTc shown in Das et al. (2008) is not evident 
in our overall sample when using the ctc — P relation (Seigar et al. 2008). However, the trend is seen when those galaxies with clearly 
fast bars are isolated. 
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Figure 10. A plot of bar strength, Qs, versus de Vaucouleurs morphological type index for a sample of four sets of galaxies including 
those listed in: (1) Table 1 and their corresponding Qg, which may be treated as an upper limit estimate of Qb\ (2) Table 1 of Rautiainen 
et al. (2008; sans those included in Table 1 of this paper) and their corresponding Qb (Buta et al. 2005); (3) Table 2 of Rautiainen et al. 
(2008) where Qb was estimated by Buta et al. (2005); and (4) Treuthardt et al. (2008, 2009) where Qg was estimated by the authors. 
The values of T were taken from HyperLeda for all of the galaxies. 



morphologies resembling that seen in the CGS B-band im- 
ages by primarily varying Q.p, a key dynamical component, 
or similarly TZ. The results of these morphological matches 
allowed us to discern those galaxies with convincingly fast 
rotating bars and correspondingly low central density dark 
matter halos. 



Our results indicate that a wide range of Mbh exists 
in both fast and slow bar galaxies. This appears to support 
the theoretical work of Booth & Schaye (2010) where the 
authors claim that the scatter in the empirical relationship 
between SMBH and halo mass (Bandara et al. 2009) is due 
to the dark matter halo concentration. Additionally, we find 
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Figure 11. Plot oiTZ = RcrI Rbar versus T with the horizontal dashed line corresponding to 7^ = 1.4, or the boundary between the fast 
and slow bar domains. The points are color coded as in Figure 10 where open circles and triangles correspond to galaxies with slow bars, 
filled circles and triangles to fast bars, and filled squares and inverted triangles to intermediate bars. The circles and squares correspond 
to data from this paper. The triangles correspond to TZ data from Rautiainen et al. (2008). The values of T for all of the data points are 
taken from HyperLeda. 



that galaxies with low central dark halo densities appear to 
follow more predictable trends in P versus T and Q b versus 
T than barred galaxies in general. The empirical MBH~Mtot 
relationship (Bandara et al. 2009) also provided us with an 
estimate of Mtot for our galaxies. For a given halo mass, 
we expect the SMBHs in gala:xies with fast bars to have 
a less-than-average mass due to the low central dark halo 
density and the resulting anemic growth predicted by Booth 
& Schaye (2010). This means that the total luminous and 
dark mass of these galaxies is expected to be the minimum 
of what would be measured observationally. This claim can 
easily be tested by obtaining detailed kinematic data and 
this will be done in a forthcoming paper. 

It should be noted that the argument by Debattista and 
Sellwood (2000) correlating galaxies harboring fast bars with 
low central density dark matter halos is theoretical and has 
not yet been tested observationally due to the inherent dif- 
ficulty in doing so. Additionally, while dynamical friction 
from higher central density dark matter halos may not yet 
have slowed recently formed fast bars, it is unlikely that 
any bars in our sample are newly formed. We presume this 
because observational evidence suggests that bars are long- 
lived structures (Jogee et al. 2004; Sheth et al. 2003) and 
the galaxies in our sample are local and therefore old. We 
also acknowledge that bars may not slow due to dynami- 
cal friction alone. Models have shown that massive bulges 
can exchange angular momentum with the bar and cause 
them to slow down (Athanassoula 2003). If we assume that 
massive bulges reside in dark matter halos of high central 
concentration (e.g. Ho 2007), then angular momentum ex- 
change with both the bulge and halo would act to slow the 
bar. In this study, we are only interested in those galaxies 
with clearly fast bars where any contribution to bar slow 



down by dynamical friction via a dark matter halo is mini- 
mal. 
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